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Compressed correlation functions and fast aging dynamics in metallic glasses 
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We present x-ray photon correlation spectroscopy measurements of the atomic dynamics in a ZrgjNiaa metaUic 
glass, well below its glass transition temperature. We find that the decay of the density fluctuations can be 
well described by compressed, thus faster than exponential, correlation functions which can be modeled by the 
well-known Kohlrausch- Williams- Watts function with a shape exponent /3 larger than one. This parameter 
is furthermore found to be independent of both waiting time and wave-vector, leading to the possibility to 
rescale all the correlation functions to a single master curve. The dynamics in the glassy state is additionally 
characterized by different aging regimes which persist in the deep glassy state. These features seem to be 
universal in metallic glasses and suggest a non diffusive nature of the dynamics. This universality is supported 
by the possibility of describing the fast increase of the structural relaxation time with waiting time using a 
unique model function, independently of the microscopic details of the system. 

PACS numbers: 64.70.pe,65.60.-t-a,64.70.pv 



O 



O 
o 



> 

O 

O 

m 



X 



I. INTRODUCTION 

Understanding physical aging in glasses is essential 
at a technological level, for controlling the material 
properties and their temporal evolution, as well as at 
a fundamental level, since glasses are often considered 
as archetypes for systems far from thermodynamical 
equilibriuniii2,. However, and despite the strong efforts 
done in the field, a complete picture of the dynamics in 
the glassy state is still missing. The main difficulty arises 
from the complexity of the glassy properties. When a 
supercooled liquid is cooled below its glass transition 
temperature, Tg, the system falls out of equilibrium in a 
metastable state which depends on the previous thermal 
history, and which evolves, with waiting time, through 
different states toward the corresponding supercooled 
equilibrium liquid phase. The existence of a multitude of 
different metastable glassy states which can be explored 
through different annealing procedures or cooling rates, 
makes difficult the development of a universal theory 
for glasses or even a plain comparison between different 
experimental results. 

Physical aging in polymeric, metallic or more conven- 
tional inorganic glasses, is commonly studied by looking 
at the temporal evolution of a given observable, keeping 
fixed other external parameters, such as the temperature 
and the pressure of the system. Macroscopic quantities 
like viscosity, refractive index, volume or elastic con- 
stants, are usually the subject of this kind of studies^"—. 
In these works, the investigated quantities evolve with 
waiting time, t^, in a way which depends on the previous 
thermal history and which can be accurately described 
by stretched exponential functions^. Similar results 
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have been reported also in studies of aging effect on 
the intensity of the high frequency tail of the dielectric 
spectrEfciir— . All these works provide information on the 
equilibration toward the supercooled liquid phase, but 
they do not give any microscopic details on the aging dy- 
namics in the glassy state. This information can instead 
be obtained by the investigation, on a microscopic scale, 
of the temporal evolution of the structural relaxation 
time, T, which represents the time necessary for the 
system to rearrange its structure toward a more stable 
configuration. Its value can be measured, for instance, 
by monitoring the long time decay of the intermediate 
scattering function, f{q,t), on a length scale 27r/q, 
with the wave- vector q corresponding to the mean inter 
particle distance. The intermediate scattering function 
is related to the density-density correlation function 
and contains all the information on the relaxation 
dynamics in the system on a length scale determined by 
q. While the dynamics in the supercooled liquid phase 
has been widely investigated in the past, the difficulty of 
probing the particle-level dynamics in realistic molecular 
glasses both with experiments and simulations has 
instead strongly delayed progress in the field. In fact 
only very few studies exist of the behavior of r below 
Tjpii^ii^r— . Numerical simulations on a Lennard Jones 
glass show that for small quenching from an equilibrium 
configuration, t slows down linearly or sub linearly with 
waiting time^^ in agreement with previous results on 
polymeric glasses^ii^iii. In addition, the long time decay 
of the correlation functions can be rescaled into a single 
master curve, leading to the validity of a time- waiting 
time superposition principle in the glassy state^^i^^. 
Recently, the aging of a Mg65Cu25Yio metallic glass has 
been studied in both the supercooled liquid phase and 
the glassy state by looking directly at the evolution of 
the dynamics at the interatomic length scale^^. The 
glass transition has been found to be accompanied by 
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i) a crossover of the intermediate scattering function 
from the well characterized stretched exponential shape 
- thus described by a stretching exponent /3 < 1 - in 
the supercooled liquid phase, to an unusual compressed 
exponential form - defined by a parameter /? > 1 - in the 
glassy state and ii) to the presence of two very distinct 
aging regimes below Tg. In particular, a fast, exponential 
growth of r has been reported for short waiting times, 
followed by an extremely slow aging regime at larger ■ 
While the latter slow regime has been already reported 
both in numerical simulations^^ and in experiments on 
polymeric glasses^^i^iii, the compressed shape of the 
correlation functions and the fast aging regime had not 
been reported before in structural glasses. Two similar 
aging regimes have been recently observed in a numerical 
study of a concentrated system of nearly hard spheresSS. 
In this case, however, the correlation functions display a 
more stretched exponential behavior, and no signature 
of a compressed shape was reported^*^. 
Compressed correlation functions and different aging 
regimes have instead been reported for out of equilibrium 
soft materials, where the dynamics is controlled more 
by ballistic-like rather than by diffusive motions2Si21. 
The previous work on the metallic glass thus suggests a 
common nature of the dynamics in metallic glasses and 
soft materials. A very similar faster than exponential 
decay and strong aging behavior for the intermediate 
scattering function has been observed also in a Zr-based 
metallic glass by Leitner and coworkers^'*. In this case, 
however, this unusual dynamics has been associated 
with the onset of a crystallization process. There are 
consequently several questions to clarify: i) are the 
results found in Mg65Cu25Yio really a universal feature 
of metallic glasses or do they depend on the microscopic 
details of the system? ii) is it possible to find a common 
description of physical aging for different glasses? iii) 
what is the effect played by the wave- vector, q, on the 
dynamics? 

Here, we answer to these questions by presenting a 
detailed investigation on physical aging in a Zr67Ni33 
metallic glass at the atomic length scale. We find 
that in this case the dynamics in the glassy state is 
also characterized by compressed correlation functions 
and different aging regimes, supporting the idea of 
a universal mechanism for physical aging in metallic 
glasses. This hypothesis is moreover strengthened by 
the possibility of describing the aging dynamics using 
the same phenomenological exponential growth and 
rescaling the results obtained for different systems into 
a single master curve. 



II. EXPERIMENTAL DETAILS 

The ZrerNiaa samples {Tg — 647 K^S) were produced 
as thin ribbons by melt-spinning at the Polytechnic Uni- 
versity of Catalonia. The Zr and Ni pure metals were 



pre-alloyed by arc-melting under Ti-gettered Ar atmo- 
sphere. The melt was then fast-quenched with a cooling 
rate of 10^ K/s by injecting it on a Cu wheel spinning 
with a 40m/s perimeter velocity. The resulting metallic 
ribbons have a thickness of 22±4 /im, close to the optimal 
value required to maximize the scattered intensity in the 
experiments. The ribbons were then held in a resistively 
heated furnace covering a temperature range between 270 
K and 1000 K with a thermal stability of less than 1 K. 
The temperature of the sample was held at T = 373 K, 
reached by heating from room temperature with a fixed 
rate of 3 K/ min. 

X-ray photon correlation spectroscopy (XPCS) measure- 
ments were carried out at the beamline ID 10 at the 
European Synchrotron Radiation Facility in Grenoble, 
France. An incident x-ray beam of 8.0 keV (A = 1.55 
A) was selected by using a single bounce Si(l,l,l) crys- 
tal monochromator with a bandwidth of AA/A « 10"^. 
Hard X-rays originating from higher order monochro- 
mator reflections were suppressed by a Si mirror placed 
downstream of the monochromator. A partially coher- 
ent beam with ~ 10^° photons/s was then obtained by 
roUerblade slits opened to 10 x 10 /im, placed ^^0.18 m up- 
stream of the sample. In order to enhance the scattered 
intensity, speckle patterns were collected in transmission 
geometry by two IkonM charge-coupled devices from An- 
dor Technology (1024 x 1024 pixels, 13 x 13 fj,m pixel 
size each one) installed in the vertical scattering plane 
^70 cm downstream of the sample. All pixels of the two 
CCD detectors were considered to belong to the same q 
with a resolution of Aq — 0.04 A~^. The detectors were 
mounted on a diffractometer and could rotate around the 
center of the sample position, covering a q range up to 
several A~^. XPCS measurements were performed for 
different wave- vectors at and around the first diffraction 
peak of the static structure factor which corresponds to 
qp — 2.55 A""'^, thereby probing directly the temporal re- 
laxation of density fluctuations at the inter-particle dis- 
tance 2'K/qp ~ 2 A. Time series of up to ^5000 images 
were taken with 7 s exposure time per frame and were 
analyzed following the procedure described in Ref.^^. No 
beam damage was observed and the amorphous structure 
of the sample was checked during the whole experiment 
by measuring the corresponding static structure factor of 
the system. 



III. RESULTS AND DISCUSSION 

As previously discussed, physical aging implies a re- 
laxation dynamics which is not the same at all times, 
but strongly depends on the sample age or waiting time. 
The temporal dependence of the dynamics can be entirely 
captured by XPCS through the determination of the two- 
time correlation function G{ti,t2) (TTCF) which repre- 
sents the instantaneous correlation of the intensity / at 
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FIG. 1. Two-time correlation functions measured with XPCS 
in ZrerNiaa at T=373 K. Data acquisition started at to=1400 s 
(top panel) and to=49000 s (bottom panel) from temperature 
equilibration. The age of the sample at any point in the Figure 
can be taken as iage ~ to + {ti +t2)/2 s. The huge slow down 
of the dynamics at the different sample ages is well illustrated 
by the clear broadening of the intensity around the diagonal 
from the left bottom corner to the top right corner in the 
bottom panel with respect to the case corresponding to a 
much shorter (top panel). 



two times ti and t2, being: 



{mm). 
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Here, the average is performed on all the pixels of the 
CCD, that correspond to the same wave-vector (^^—. 
Figure [1] shows the TTCF measured in ZrgyNiaa in the 
glassy state at two different sample age times from tem- 
perature equilibration. The colors represent the values 
of G. The diagonal from bottom left to top right cor- 
responds to the time line of the experiment and it thus 



indicates the sample age. The width of this diagonal 
contour is instead proportional to the relaxation time r. 
The steady slowing down of the dynamics, and thus the 
increase in the relaxation time due to physical aging, is 
illustrated by the broadening of the diagonal contour on 
increasing the sample age. This effect is clearly visible 
in the top panel of Figure [T] where we report data taken 
between ~ 1000 and ^ 5000 s after temperature equili- 
bration. After almost 13 hours of annealing (lower panel 
of Figure [J), the dynamics was so slow that the inten- 
sity along the diagonal covered almost the whole image, 
despite the time interval was almost the double with re- 
spect to the one reported in the top panel. After this 
very long annealing the width of the diagonal contour 
seems to be constant with time, as it would be in the 
case of stationary equilibrium dynamics. As discussed 
later, here the effect of physical aging is probably so slow 
that the dynamics appears as stationary at least on our 
experimental time scale. 

Although the dynamics is not stationary, the aging is suf- 
ficiently slow that it is possible to find time intervals over 
which the G{ti,t2) can be averaged and thus to get the 
standard intensity autocorrelation function 
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where (...) is the temporal average over all the possible 
starting times ti . The intensity autocorrelation function 
is related to the intermediate scattering function through 
the Siegert relation g2iq,t) = 1 + B{q) \f{q,t)\^ , with 
B{q) a setup-dependent parameter—, and its determina- 
tion allows capturing snapshots of the dynamics at dif- 
ferent waiting times. The Siegert relation requires the 
temporal average to be equal to the ensemble average. 
This condition is fulfilled for systems at thermodynamic 
equilibrium. The extension to non-ergodic systems, such 
as glasses, is possible thanks to the use of an area detec- 
tor, like the CDD used during our experiments. In this 
case in fact, the intensity distribution is sampled cor- 
rectly by the average over all the pixels, independently 
of the nature of the dynamics^^i^. 

The time interval for temporal averaging was chosen 
small enough that the shape of (72(9, t) was not noticeably 
affected by the evolution of the dynamics, while keeping 
a good signal to noise ratio. For each set of images, the 
sample age was defined SLSt^ = to + (tf — ti)/2, where to is 
the delay of the starting point for the measurement from 
temperature equilibration, and t f and ti are the final and 
starting times of the interval used for the average. 
The a relaxation process in glass-formers can be well de- 
scribed by using the empirical Kohlrausch Williams Watt 
(KWW) model function^ 



f{q,t) = f,{T)-eM-{t/rr] 



(3) 



In this expression, r is the characteristic relaxation time, 
/3 is the stretching parameter which quantifies the de- 
viation from a simple exponential behavior, and fq{T) 
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FIG. 2. Selection of the intensity correlation functions mea- 
sured with XPCS at T = 373 K for different waiting times 
from temperature equilibration. On increasing the waiting 
time, the decay clearly shifts toward longer time scales due 
to physical aging. The data are reported together with the 
resulting fit obtained using equation Q. 

is the nonergodicity factor. It represents the height of 
the intermediate-time plateau before the final decay as- 
sociated to the a-relaxation and it slightly depends on 
temperature in the glassy state. 

Following the Siegert relation, the intensity correlation 
functions measured with XPCS can then be described by 

52(g, t) = l + c{q, T) ■ exp[-2(Vr)^] (4) 

where c = B{q)f^, being B{q) the parameter entering 
the Siegert relation. 

A selection of normalized intensity autocorrelation func- 
tions obtained at T = 373 K for a fixed q = Qp and differ- 
ent waiting times is reported in Figure [H The data are 
shown together with the best-fit line shapes obtained by 
using equation The intensity correlation functions 
decorrelate completely on the experimental time scale, 
implying that the system is able to rearrange its struc- 
ture on the probed spatial length scale of few A. The 
effect of aging is illustrated by the shift of the decay to- 
ward longer time scales upon increasing the sample age, 
in agreement with the analysis of the two times correla- 
tion functions reported in Figure [TJ As shown in Figure 
[21 aging is so important that the decay time increases of 
almost a factor 10 in less than 13 hours. This fast aging 
is likely to be the reason of the well known rapid embrit- 
tlement of metallic glasses^^*^. 

Interestingly, the correlation functions display a com- 
pressed, thus faster than exponential behavior, which can 
be well described by equation Q with a shape param- 
eter /3 ^ 1.8, independent of the sample age. Such a 
high value is in agreement with the results reported for 
other metallic glasses^^i^l, while it cannot be explained 
within the current theories for the dynamics in the glassy 
state'^*^"'^^. Based on the observations of the dynamics in 



the liquid above the glass transition temperature, these 
theories predict a stretched, thus with /? < 1, decay of 
the density fluctuations even in the glassy state. The idea 
is that the dynamic heterogeneities present in the liquid 
phase partially freeze below Tg, leading to a narrower 
distribution of relaxation times, and thus to an increase 
of the shape exponent, up to the limiting case of /3 = 1. 
This explanation well describes the correlation functions 
calculated in numerical simulations^ and the increase of 
(3 measured in some polymeric glasseai^iii, but cannot 
take into account the behavior reported in Figure [H 
From the comparison with other systems, it seems then 
clear that metallic glasses behave in a different universal 
way, independently of their atomic structure. In the case 
of Mgg5Cu25Yio, the unusual faster than exponential be- 
havior has been attributed to the strain field generated by 
a random distribution of slowly-evolving sources of inter- 
nal stresses stored in the system during the quench, and 
which can then be released upon annealing the system 
close to Tg"^^-^^. This interpretation suggests the exis- 
tence of a different kind of dynamics, rather than pure 
diffusion, in metallic glasses, which display many similar- 
ities with that of out-of-equilibrium soft materials, such 
as concentrated colloidal suspensions and nanoparticle 
probes in a glass former matri x^^i'^^" — . 
In the case of ZrerNiaa the shape parameter is even larger 
than the one reported for Mg65Cu25Yio (/3 ~ 1.5 for 
a similar thermal path). This may suggest the pres- 
ence of larger atomic mobility or more internal stresses 
in ZrgyNiaa , a conclusion which however requires further 
investigations. While the value /3 ^ 1.5 has been often 
observed in numerous experimental studies of soft mate- 
rials and corresponds to the early-time regime proposed 
in the model of Bouchaud and Pitard^-, values of (3 as 
large as those here reported for ZrerNisa have been ob- 
served only in another Zr-based metallic glass^ and in 
few other soft systems, and their interpretation is still far 
from being reached^i^. 

The waiting time dependence of the structural relax- 
ation times obtained from the analysis of the g2{%,i) is 
shown in Figure [3] a) (squares) . Data corresponding to 
different q values around the main peak of the diffraction 
pattern are also reported. As explained in the experimen- 
tal details, the low scattered signal allowed us to measure 
just one q vector for each detector setting. Consequently 
the different qs are taken at different tw from temperature 
equilibration. At all the investigated qs, r grows rapidly 
with the sample age, in agreement with the results found 
in the Mg65Cu25Yio glasses.. Interestingly, the dynam- 
ics is characterized by structural rearrangements which 
require ~ 10'^ — 10^ s thus, relatively short times with 
respect to the common idea of an almost frozen dynam- 
ics in the deep glassy state. These values are however in 
perfect agreement with those reported for both the Mg- 
based and Zr-based metallic glasse o^^'^*^ and suggest the 
presence of a peculiar dynamical behavior at the atomic 
length scale. Hence, it would be interesting to investigate 
the nature of the atomic dynamics in the glassy state 
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FIG. 3. a) Waiting time dependence of the structural relax- 
ation time measured at r=373 K for different wave-vectors 
around the first sharp diffraction peak (g^ = 2.55 A^^). The 
inset b) reports a zoom of the data taken at large waiting 
times, c) Corresponding shape parameter, d) Selection of 
the correlation curves corresponding to the data reported in 
a) as a function of the time rescaled by the structural relax- 
ation time obtained from the analysis of the curves with the 
KWW model function of Eq. @. The inset e) indicates the 
position of the different investigated qs. 

through the q dependence of the structural relaxation 
time. Unfortunately, this information would require at 
each investigated q a good knowledge of the aging in a 
wide range, which can be gained only by measuring 
the dynamics simultaneously for different wave-vectors. 
A higher incident flux, as the one available at free elec- 
tron laser sources, will make this kind of studies possible. 
Interestingly, the shape parameter does not display any 
dependence on q ot tw, and remains constant with /3 ^ 
1.8 ± 0.08 for all the investigated qs (Figure E] c)). This 
means that it is possible to rescale all the data sets to a 



single master curve and define a time- waiting time- wave 
vector superposition principle in the glass. This scaling 
is reported in |3] d) where a selection of intensity auto 
correlation functions corresponding to the data reported 
in panel a) are plotted as a function of t/r. In the case 
of Mg65Cu25Yio, /3 was found to be temperature and 
independent in the fast aging regime^^. Our results indi- 
cate that this parameter is also not affected by the choice 
of the wave- vector, at least for values of q close to that of 
the maximum of the static structure factor (see Figure [3] 
e)). 

As shown in Figure [3] a) and b), for q = qp and large 
t^u, the fast aging regime seems to change to a more 
stationary or very slow one, as also suggested by the 
analysis of the two times correlation function reported 
in the bottom panel of Figure [1] The presence of a com- 
plex hierarchy of aging regimes has been observed also in 
Mg65Cu25Yio but for temperatures very close to the glass 
transition one {T/Tg = 0.995). The results here reported 
for T/Tg = 0.577 show that this second aging regime can 
be found also in the deep glassy state, after a sufficiently 
long anneahng time. In the case of Mg65Cu25Yio, the 
second aging regime was accompanied by a decrease in 
the shape parameter. As shown in Figure |3]c), in the 
present case the lack of a full complete decay of g2 (t) for 
the highest leads to a larger error in the determination 
of /3, which can be determined only for > 50000 s and 
remains however similar to the value obtained at lower 

From the energetic point of view, the abrupt crossover 
to an almost stationary dynamics suggests the exis- 
tence of a deep energy minimum in the potential en- 
ergy landscape, where the system remains trapped for 
a very long time before eventually jumping into a dif- 
ferent configuration^i^. A similar plateau has been re- 
cently observed in the enthalpy recovery of glassy poly- 
mers and has been interpreted as a thermodynamically 
stable stateiS. This explanation cannot be extended to 
the case of metallic glasses. In fact, the comparison be- 
tween XPCS and low temperature equilibrium viscosity 
data for Mg65Cu25Yio suggests instead a gradual ap- 
proach toward equilibrium^^, in agreement with the sub 
linear aging regime found in other polymeric glasses and 
coUoids^i^ii^. 

Following the procedure reported in Ref<^, we can de- 
scribe the short tw behavior of Tq^{tw) using the phe- 
nomenological expression 



'(r,i^) = To(r)exp(t„/T*) 



(5) 



where tq{T) is the value of r for = at the temper- 
ature equilibration time, and r* is a fitting parameter 
which describes the rate growth and is independent of 
temperatureSi. A fast exponential growth with tyj is ev- 
ident also at the other investigated qs. However, in that 
cases the lack of information on the behavior at shorter 
tu, does not allow for a more quantitative description of 
the aging dependence for r. 
For q = qpwe find tq = 407 ± 6 s and r* = 6490 ± 170 
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FIG. 4. Structural relaxation times measured in two as 
quenched glasses of Mg65Cu25Yio at T/Tg = 0.995 with dif- 
ferent thermal histories and in ZrerNisa at T/Tg = 0.577. All 
the data are normalized for the tq parameter obtained from 
the analysis of the fast aging with the empirical model func- 
tion ((5| and are reported as a function of the time rescaled 
by the r* parameter. The data for Mg65Cu25Yio are taken 
from ReiM. 



s. The latter value is quite close to the one reported 
in leiM for MgesCuasYm at T/Tg ^ 0.995 and it is 
of the same order of magnitude of that previously re- 
ported in the case of colloidal gels and Laponit o^^'^^ , thus 
strengthening the idea of a universal out-of-equilibrium 
behavior. The comparison with the Mg-based metallic 
glass suggests that the physical aging is very similar for 
both systems and that it should be possible to rescale 
the rapid growth of r(t^) to a unique curve, indepen- 
dently on the microscopic details of the system. This 
scaling is shown in Figure |31 where the structural relax- 
ation times of ZreyNisa at T/Tg = 0.577 are compared 
to those measured in Mg65Cu25Yio at T/Tg = 0.995 and 
for two different thermal paths. While squares and circles 
correspond to as-quenched samples heated up to a given 
temperature, stars refer to a state reached after having 
partially annealed the glass at higher temperatures with- 
out equilibration^^. All the data are rescaled for the cor- 
responding To parameter obtained from the analysis of 
the fast aging using equation ([SJ and are reported as a 
function of the time rescaled by the r* parameter. For 
short t^, all the data fall into a unique master curve, 
suggesting a common origin of the fast aging regime in 
metallic glasses. Differently, the crossover to the subse- 
quent almost stationary dynamics strongly depends on 
the temperature and thermal path, and is reached only 
for very long annealing times (circles) or after particular 
thermal treatments (stars). 



IV. CONCLUSIONS 

In conclusion, we have studied the physical aging in 
a ZrerNiaa metallic glass, by following the temporal 
evolution of the structural relaxation process at the 
interatomic length scale. We find that the decay 
of the density fiuctuations can be well described by 
compressed, thus faster than exponential correlation 
functions, suggesting a peculiar origin for the atomic 
dynamics^i^. This behavior seems to be universal 
in metallic glasses^^iM, and cannot be explained by 
the well-known theories for glasses dynamics, which 
predict a stretched exponential shape of the correlation 
functions, even in the glassy stated"—. Model systems 
such as hard spheres and Lennard- Jones are unable to 
reproduce these intriguing dynamic effectai^i^S. In these 
model systems, indeed, the dynamics appears to be 
always characterized by a broad spectrum of relaxations. 
Conversely, very similar results have been previously 
reported for out of equilibrium soft materials, suggesting 
a universal microscopic dynamics for these very different 
systemsSS,. In these systems, the compressed behavior 
arises from ballistic-like motion due to the presence 
of internal stresses in the out of equilibrium stated!. 
It would be interesting to check whether this is also 
the case for metallic glasses. However, a q dependence 
study would require a much stronger signal than the one 
achievable even at a third generation synchrotron source 
as the one used for the present study. Notwithstanding, 
indication of the presence of distinct processes of atomic 
motion in metallic glasses have been reported for differ- 
ent Zr-based metallic glass-formers studied with nuclear 
magnetic resonance^. In these systems the diffusive, 
collective motion is found to be accompanied by hopping 
which becomes dominant in the glassy state, and could 
then be related to the observed compressed relaxation 
process in our glass. 

Albeit we could not get information on the q dependence 
of the relaxation time due to the impossibility of mea- 
suring simultaneously different gs, we have investigated 
the nature of the dynamics at different wave- vectors and 
sample ages. We find that the shape of the correlation 
functions does not show any dependence neither on 
the sample age nor on the wave- vector, at least for qs 
around the maximum of the static structure factor of the 
system, which is at qp = 2.55 A"-'^. These findings allow 
then to rescale all the correlation curves on the top of 
each other and define a time-waiting time-wave-vector 
superposition principle in the fast aging regime. 
The value of /3 seems to depend on the details of the 
system, being larger in Zr-based systems than in the 
previously reported Mg65Cu25Yio^'^. As both these 
glasses were obtained with the same thermal protocol, 
we speculate that the observed difference could be re- 
lated to a larger atomic mobility in the first system and 
to the consequent presence of larger internal stresses. 
Despite this small difference, the comparison between 
the two metallic glasses suggests the existence of a 
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universal dynamical behavior, also indicated by the 
scaling of the fast aging regime for the relaxation time to 
a single curve, independently of the microscopic details 
of the systems. Additional theoretical and experimental 
works is clearly required to elucidate this remarkable 
dynamical behavior of metallic glasses. 
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